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Abstract. An effective field theory developed for systems interacting through 
CJ ■ short-range interactions can be applied to systems of cold atoms with a large 

^ . scattering length and to nucleons at low energies. It is therefore the ideal 

tool to analyze the universal properties associated with the Efimov effect in 
three- and four-body systems. In this progress report, we will discuss recent 
^ . results obtained within this framework and report on progress regarding the 

! inclusion of higher order corrections associated with the finite range of the 

underlying interaction. 
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1 Introduction 

Universal properties are one reason why few-body systems with a large scattering 
length a have received a lot of attention lately. For example, a two-body system 
with a large positive scattering length a has a bound state with binding energy 



where m denotes the mass of the particles (note that we set h = 1 for con- 
venience). In the three-body sector a large two-body scattering length leads to 
further universal features known widely as Efimov physics [1]. These are fasci- 
nating and exciting since they relate physical effects that occur at very different 
length scales. In experiments with ultracold atoms Feshbach resonances allow 
the interatom scattering length to be tuned from — oo to oo, while also the in- 
ternucleon interaction exhibits a scattering length large compared to the range 
of the interaction. 

An understanding of such systems improves therefore also our understanding 
of low-energy nuclear physics. Moreover, a simple approach only in terms of the 
large scattering and the other remaining low-energy parameters of the effective 
range expansion grants also the possibility to address important problems in 
nuclear few-body systems in a very effective manner. One prominent example is 
the model-independent description of electroweak reactions relevant to nuclear 
astrophysics. An approach whose parameters are only related to the effective 
range expansion might seem inappropriate for such nuclear reactions but it is 
important to realize that many important processes occur at energies well below 
the energy scale set by the range of the interaction (i.e. the pion mass m,r)- 

Since such an approach will fail for observables for which the pion-exchange 
character of the internucleon interaction is important, it can therefore also answer 
the question for which measurable quantities the chiral character of the nucleon- 
nucleon interaction is relevant. This addresses the question what low-energy in- 
formation in particular is required to describe nuclear systems at different energy 
scales. 

The framework that we will use to address processes in such systems is that 
of effective field theory (EFT). It allows to obtain model-independent results 
in a small parameter expansion. The advantages of a calculation of observables 
in the EFT framework are numerous. Firstly, it allows for a clear separation 
of the unknown ultraviolet and known low-energy physics which eliminates any 
possible model-dependence. Secondly, the accuracy of calculated observables can 
be systematically improved by calculation of another order in the low-energy 
expansion. And thirdly, the quantum field-theoretic setup brings along the ad- 
vantages of mathematically tools successfully used in high-energy physics, such 
as Feynman diagrams, renormalization and regularization. 

Here we will concentrate on applications of one particular EFT (that we will 
call short-range EFT) that is tailored to the problem of short-range interactions. 
The expansion parameters in this EFT are R/a and kR, where R denotes the 
range of the underlying interaction and k the momentum scale of the process un- 
der consideration. At leading order, this EFT framework allows us to analyze the 
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properties of systems interacting through zero-range interactions. Finite range 
effects can then be accounted for by calculating subleading orders in the EFT 
expansion. 

The short-range EFT provides therefore a good framework to analyze the 
few-body dynamics in resonantly interacting ultracold gases. Binding energies 
and scattering observables can be calculated and the experimental signals of 
Efimov physics can therefore be identified. 

Predictions for low-energy scattering processes and binding energies can also 
be made in the nuclear case. However, the field-theoretic formulation of this 
framework facilitates in addition a calculation of electroweak observables relevant 
to nuclear astrophysics as discussed above. 

The purpose of this review is to display the similarities between low-energy 
atomic and nuclear physics and to demonstrate that the same framework can 
be applied to both sectors to obtain relevant and interesting results. In con- 
tradistinction to recent reviews on similar topics [21 El HI [5] , we will present the 
applications and results for atomic and nuclear systems alongside. By doing this 
we want to emphasize the similarities between low-energy nuclear physics and 
the few-body dynamics in resonantly interacting ultracold gases of alkali atoms. 

We will be mainly interested in few-body systems and will therefore mainly 
focus on systems that can be described either as two, three or four-body systems. 
It should also be noted that this is not a complete review of work done for systems 
with short-range interactions but instead it should be considered as a progress 
report that gives a brief introduction to the short-range EFT, that presents recent 
calculations performed in the few-body sector and discusses then an incomplete 
list of problems that should be addressed in the near future. 

In the following chapter, we will thus give a brief introduction to the concept 
of EFTs and will discuss in detail the short-range EFT whose application to 
different problems will be covered in the remaining sections of this review. In 
section [3l we will give a review of recent applications of the short-range EFT 
to systems of ultracold atoms. In section [J] we will review results obtained with 
short-range EFT for nuclear systems such as the three- and four-nucleon systems 
but also heavier systems such as halo nuclei. We will conclude with a summary. 

2 Effective Field Theories 

Effective field theories are an excellent tool to describe the interactions in systems 
with a separation of scales. The usual example of such a scale separation is the 
ratio of a small momentum over a heavy mass scale. The weak interaction for 
example, can be described without explicit use of heavy vector bosons if the 
energies of the processes are well below the mass scale of the exchange particles. 
In this case the weak interactions can be described very well with the familiar 
Fermi interaction. 

Here, we will focus on the short-range EFT, which is tailored to describe 
non-relativistic particles with a large scattering length a at low energies. We will 
discuss the Lagrangian which is at the heart of the EFT and describe how it can 
be used along with an ordering scheme called power counting and Feynman rules 
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to obtain a powerful framework for the computation of low-energy observables 
of few-body systems with a large two-body scattering length. 

Very good introductions to the general concepts of EFTs have been written 
by Kaplan [6], Lepage [7| and Polchinski [8j. An excellent introduction to EFTs 
applied to nuclear systems has been given by Phillips [9]. 



2.1 The Effective Range Expansion 

We will consider identical particles (bosons for simplicity) interacting via a po- 
tential that has a finite range R. In this case the amplitude for the scattering of 
two particles can be written as 



An SvTv-^ 21 + 1 ^. 

A{k,cosO) = — > r — ii(cosi 

m ^-j^ k cot oi — ik 



where k denotes the relative momentum between the two particles, m is the mass 
of each of these particles and (5/ denotes the scattering phase shift in partial wave 
/. 

At low momenta, the expression in the denominator can be expanded in even 
powers of the momentum k 

k^i+^ cot 5i = -- + -k^ + Oik"^) . (3) 
ai 2 

For the S- and P-wave phase shifts, the first terms in the expansions are 

11 3 

k cot 6o = 1- -rgk'^ + ... , k cot 6i = --ttt^ + ... , (4) 

a 2 k'^a^ 

where a, ap are the S- and P-wave scattering lengths, respectively, is the 5- 
wave effective range and the ellipses denote higher order terms in the expansion. 

For large positive scattering length, one finds a bound state as discussed 
above and we will frequently use the effective range expansion (ERE) around 
the bound state pole: 

/c cot (5o = -7 + ^(7' + fc') + ... , (5) 

where Ed = 7^/m and p = up to the orders considered in this work. 

The low-energy physics encoded in the effective range parameters depends on 
the relative size of these parameters. In the following we will consider two special 
cases which we will refer to as the natural case and the unnatural case. In the 
natural case, all effective range parameters are of natural size, i.e. of the order 
of the range of the interaction a ~ ~ i?. In this case the scattering amplitude 
in Eq. ([2]) can be expanded in powers of k: 

j^^k) = -^\l-iask + iasrs/2-a'',)k'' + ...] . (6) 
m 

For simplicity we have considered here only the 5- wave projected part of the 
scattering amplitude. This expansion will converge for momenta <C 1/R . 



6 



Low-Enorgy Universality in Atomic and Nuclear Physics 



Now consider the case in which the scattering length a is much larger than 
all other effective range parameters \a\ ^ \rs\ ~ R. This separation of scales 
indicates the presence of non-perturbative physics through a (virtual) bound 
state with energy ^ l/(mo^). The expansion above is only valid for ak <C 1 
which reduces the radius of convergence dramatically and renders it useless in 
the limit of infinite scattering length. We will instead expand in powers of kR 
and keep ak to all orders: 



^n^ ^TT 1 
A{k) ■■ 



m 1/a + ik 



rJ2 ,2 irs'2f ,4 
la + ik {la + ikY 



(7) 



The resulting expansion still converges for kR <^ 1 and reflects the fact that 
systems with a large two-body scattering length require a non-perturbative re- 
summation at leading order. 

2.2 The Short-Range EFT 

Our goal is to describe the dynamics of particles interacting through a short- 
range potential with an EFT. The Lagrange density for an EFT is generated 
by writing down all possible operators built from the available degrees of free- 
dom in accordance with the required symmetries. At sufficiently low momenta, 
non-relativistic particles can be described by an EFT built up from contact inter- 
actions alone. This EFT can be applied if the relative momenta k of the particles 
is much smaller than the inverse of the range of the underlying interaction R: 

k<^l/R. (8) 

The convergence radius of our EFT will therefore agree with that of the ERE 
and we can therefore anticipate that we will be able fix any free parameters in 
the two-body sector from the ERE. 

The only degrees of freedom we require for our EFT are therefore the atoms 
themselves. Having identified the relevant degrees of freedom, we generate the 
effective Lagrangian by writing down all possible operators satisfying the con- 
straints of Galilean, parity and time-reversal invariance and locality: 



The ellipses represent operators of higher dimension which means terms with 
more derivatives and/or more fields. We have neglected relativistic effects which 
are suppressed by factors of (p/M)^. Dq and £^0 denote the leading three- and 
four-body interactions. 

Every EFT contains therefore an infinite number of two- and many-body 
operators. This framework might therefore seem useless since an infinite number 
of Feynman diagrams arises for every observable and infinite number of coupling 
constants needs to be determined from experiment. Here, the afore-mentioned 
separation of scales allows us to resolve this problem. The separation of scales 
usually leads to a ratio of a light scale over a heavy scale (e.g. a small momentum 
divided by a large mass) which is exploited as a small expansion parameter ckeft- 
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Figure 1. The sum of diagrams including the Co vertex contributing to two-body scattering. 

Appropriate renormalization of the vertex constants (also known as low-energy 
constants) will then lead to an ordering of all possible Feynman diagrams (power 
counting) such that only a well-defined number of low-energy constants are re- 
quired at each order of the expansion of a matrix element in powers of aEFT- 
This power counting arises in the short-range EFT by demanding that the scat- 
tering amplitude for two-body scattering reproduces the appropriate momentum 
expansion (Eqs. ([6]) and ©) of the scattering amplitude. 

Feynman rules: 

As in standard quantum field theory, amplitudes for processes are calculated 
by using Feynman rules derived from the Lagrange density. The Feynman rules 
for the short-range EFT are particularly simple 

• Assign non-relativistic four-momenta {po , p) to all lines and enforce mo- 
mentum conservation at each vertex. 

• For each vertex include a factor —i times the low-energy constant of the 
corresponding operator. 

• For each internal line include the propagator with four-momentum {po,p) 

iS{q,,qo) = -2 , (9) 

(10- 

• Integrate over all undetermined loop momentum using the measure 

J 27r (27r)3 ' 

The energy integral can be evaluated using contour integration. 

• Multiply by a symmetry factor 1/n! if the diagram is invariant under the 
permutation of n internal lines. 

Power counting 

It is clear from the discussion on the expansion of the two-body amplitude 
for large scattering length that the power counting requires a resummation of 
an infinite set of diagrams to obtain the pole structure of Eq. ([7]). Here, we will 
verify by explicit calculation that the leading order two-body amplitude can be 
reproduced by summing up all diagrams including only the Cq vertex as shown 
in Fig. [TJ 
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Using the Feynman rules given above, the sum of diagrams shown in Fig. [T] 
can be expressed as 



oo 

A{E) = -iCo^(-iCoT/2r 

n=0 

-iCo 



1 + iCoI/2 ' 

where I is given by 

dqo dq^ i 



(10) 



2^ (2^)3 _ go _ |i + ie go _ |i + ie 



This integral has been regularized with an ultraviolet cutoff A. Inserting Eq. (112 
in (jlOp . we end up with the expression 



. Svr 
A = — 

m 



Svr 2A , 

h V-n^Po 



mCo vr 

We reproduce the leading order term in Eq.Q by setting 



(13) 



ml — zAa/TT 

In the two-body system the leading order amplitude is therefore given by the sum 
of all diagrams which contain only the Cq vertex. Subleading corrections can then 
be calculated from perturbative insertions of the higher order operators dressed 
to all orders by Cq. 

The short-range EFT power counting for the two-body sector was developed 
in Refs. [ini [HI 111 HSl [H] . A regularization scheme alternative to cutoff regular- 
ization which makes the powercounting on a diagram-by-diagram basis explicit 
was developed in Refs. pTl [T2] . 

It has to be emphasized that the powercounting is of course different in the 
case of a scattering length of natural size, i.e. a ~ i?. In this case the pow- 
ercounting and therefore the order at which a particular diagram contributes 
follows directly from naive dimensional analysis of the corresponding operators. 
This EFT is only of limited interest in the few-body sector, but it has found 
interesting applications in the many-body sector [151 HB] • 



2.3 The Three-Body System 

Applications of the framework laid out above to the two-nucleon system have 
been very successful (for a review see Ref. [1]). 

Bedaque, Hammer and Van Kolck were the first to consider the three-body 
system using the short-range EFT [l7\ |18] . They used an equivalent form of the 
Lagrange density given in Eq. ([9]) that turns out to simplify the treatment of the 
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Figure 2. The dressed dimeron propagator 



three-body system. This form uses an auxihary field, the dimeron T, which has 
the quantum numbers of a two-body state |17j : 

£ = ,/,t Lq^^^] t/j-^AT^T- -^(TVV' + h.c.) + /iT'frV'V---- (15) 
y 2m J V2 

The Lagrange density above is equivalent to the density in Eq. ([9]) if the low- 
energy constants are chosen to be 2g^ j A = Cq and —IShg'^/A-^ = Dq. The 
dimeron field T has the bare propagator 

It is clear that the power counting introduced in the previous section dictates 
a resummation of the particle-particle loops as shown in Fig. [2l We therefore 
obtain the dressed dimeron propagator 0: 



4tt i 



mg^ -l/a + yJ-mpQ + p2 



(17) 



le 



This propagator has a pole at pq = p^ /{Am) — {AiTA/{mg'^))'^/m. Evaluating the 
residue of Eq. ()17p leads to the wavefunction renormalization factor 

Zd = . (18) 

Let us consider now elastic scattering between an atom and a dimer shown 
diagrammatically in Fig [3l Using the Feynman rules discussed previously with 
the added information on the dimeron propagator, we can consider the amplitude 
for elastic particle-dimer scattering. 

t(p,k;S) = 5 - + h 



(27r)3 + ^-iq^l^-mE - it 



i(q,k;£;) ^ h 



-mE -\- + p^ + p ■ q 2mg^ 



(19) 



The reason why we have included the three-body interaction h will be discussed 
below. 



^No cutoff dependence appears in the expressions for the two-body propagator in Eq. (|17p 
since we assume that the bubble loop integrals from the two-loop sector were evaluated using 
dimensional regularization with minimal subtraction. 
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Figure 3. The integral equation for the atom-dimer scattering amphtude t 



In order to relate the amplitude to observables, we have to multiply t with the 
wavefunction renormalization Z/j factor from Eq. (jlSp . It is therefore convenient 
to rewrite the integral equation: 



«ii(p,k;^) 



+87r 



167r 
ma 



1 



^2 _|_ p2 _ ^JJJ _|_ p . 
1 



+ 



HjA) 



(27r)3 + y^3q'^/4-mE -ie 



t/?(q,k;S) 



-mE + + + p • q 



+ 



HjA) 



(20) 



where tji = Zot and h = 2mg^H{A)/A^. The total energy in this process is given 
by -E = — Ef) and the on-shell point is given by k = p. It is straightfor- 

ward to decompose the amplitude into contributions from channels with different 
orbital angular momentum quantum number /: 



1 

ti{p,k;E) = -J ^dx Pi (x) t{p,k;E) . 



(21) 



where x = p • k/(p/c) and Pi{x) denotes a Legendre polynomial. Projecting onto 
S-waves {I = 0) gives 

to{p,k;E) 



Stt 
ma 



pk 

2 

+- 
TT Jo 



pk + k'^-mE\ 2H{A) 



pk + k^ 



+ 



dqq 



mE J ' yl2 
t{q,k;E) 



-l/a + i/3g2/4 - mE - ie 



1 , [ p^ +pq + q^ 
— m 
pq 



mE 



■ pq + q^ — mE 



+ 



2H{A) 



(22) 



Here a cutoff A has been introduced to make the integral equation well-defined. 
Equation ([22|) is then related to the atom-dimer phase shift via 



to{k,k) 



Sir 



1 



m k cot 6 AD — ik 



(23) 



Equation (|20|) (without the three-body force) is also known as the Skorniakov- 
Ter-Martirosian (STM) equation, named after the first ones to derive an integral 
equation interacting through zero-range two-body interactions [19j . 

Bound states can be found by solving the homogeneous version of the integral 
equation in Eq. (I22p . So far we have not explained why we included the three- 
body force in our leading order equation. Naive dimensional analysis by counting 
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Figure 4. The shallowest three-body binding energies indicated by the solid, dashed, and 
dash-dotted lines as a function of the momentum cutoff A. 

the mass dimension of the three-body operator suggests that the three-body force 
is a higher order effect. 

Let us therefore consider the bound state problem and set the three-body 
force to zero. The resulting strong cutoff-dependence of three-body binding en- 
ergies is shown in Fig. HI This cutoff-dependence is no residual dependence which 
becomes weaker with increasing cutoff but is a genuine result of the use of zero- 
range interactions in the three-body sector. It was in fact already pointed out by 
Danilov in the 1960's that the STM equation has no unique solution ^0] and it 
was later suggested by Kharchenko to fix the cutoff and to treat it as a parame- 
ter [21]. The field-theoretic perspective implies that three-body observables are 
sensitive to short-distance effects that have not been properly renormalized. 

Bedaque, Hammer and van Kolck [171 [T8] showed that the inclusion of the 
three-body force as in Eq. (j22p results in fully renormalized observable^. The 
three-body force H(A) is therefore included and its value is determined by ad- 
justing Eq. (j22|) to a given three-body datum (such as the binding energy of a 
three-body bound state). It was shown furthermore that the cutoff dependence 
of H{A) can be approximated with 



^An alternative perspective can be won by an analysis of this problem in coordinate space. 
An analysis in the hyperspherical formalism shows that the three-body problem with infinite 
scattering length becomes a Schrodinger like problem with a 1 /TZ^ potential, where 71 denotes 
the three-body hyperradius. The solution to this problem requires a short-distance boundary 
condition which is determined by a three-body datum. An excellent discussion of the analysis 
of this problem in the hyperspherical formalism is given in [22] . 



H{A) 



sin (sq ln(yl/L3) — arctan(l/so)) 



(24) 



sin (so ln(yl/L3) -|- arctan(l/so)) 
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Figure 5. Coupling constant as a function of the cutoff A. The solid shows a fit to Eq. (|24|l 
for an arbitrarily chosen three-body renormalization point. 



where -L3 is a parameter fixed from experiment and sq ~ 1.00624. It can easily 
be seen tliat Hi^A) is periodic in tlie cutoff. A rescaling of the cutoff yl by a factor 
of exp(n7r/so) ~22.7 gives back the same result for H 

H{A) = H{Aexp{mr/so)) . (25) 

This particular running of the coupling constant H is called a limit cycle and 
was first disussed as an additional type of renormalization group flow in Ref. |23j . 
The limit cycle is also reflected in observables: whenever the cutoff is increased 
by a factor of 22.7, the number of bound states in the spectrum increases by one 
(as can be seen in Fig. In the limit of infinite scattering length, one finds an 
infinite tower of three-body bound states: 

= (e-2-/so)n-n,^2/^^ (26) 

where is the binding wavenumber of the Efimov trimer labeled by n^, and 
m is the mass of the particles. The short-range EFT reproduces therefore at 
leading order the Efimov effect. Vitaly Efimov discovered in the 1970s that the 
zero-range limit of the 3-body problem for nonrelativistic particles with short- 
range interactions shows discrete scale invariance. If a = ±00, there are infinitely 
many 3-body bound states with an accumulation point at the 3-atom scattering 
threshold. These Efimov states or Efimov trimers have a geometric spectrum |Tj. 
Furthermore, he pointed out that these results were also valid for finite scattering 
length as long as a ^ . The short-range EFT is therefore an approach in which 
the properties of few-body systems related to Efimov physics can be analyzed 
but that furthermore facilitates the systematic inclusion of perturbations due to 
the finite range of the underlying interaction. 
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Figure 6. The a'^-K plane for the three-body problem. The allowed regions for three-atom 
scattering states and atom-dimer scattering states are labeled AAA and AD, respectively. The 
heavy lines labeled T are two of the infinitely many branches of Efimov states. The cross- 
hatching indicates the threshold for scattering states. 

In the discussion above, we have considered systems with identical mass 
which we will consider for the main part of this work. Nonetheless, it should 
be pointed out that the consequences of different masses on the discrete scale 
invariance/limit cycle are well- understood [2]. 

One-parameter correlations: 

The necessity of two counterterms at leading order to obtain cutoff indepen- 
dent results implies that two very different types of one-parameter correlations 
exist. One can for example consider correlations between two- and three-body ob- 
servables by changing the two-body counterterm and therefore effectively chang- 
ing the scattering length. A well-known example of such a correlation is the 
so-called Efimov plot that is shown in Fig. [6l In this plot the binding momentum 



of the nth trimer state Kn = y mE>p is plotted against the inverse scattering 
length. At the origin of the plot the three-body system shows the afore-mentioned 
exact discrete scaling symmetry with respect to the binding momenta along the 
y-axis. A change of L3 in Eq. ()24p corresponds therefore to a rescaling of the Efi- 
mov plot. Any point in this plot is therefore equally well-suited to characterize 
the short-distance behavior of the three-body amplitude and can be used in- 
stead of L3 as the three-body parameter. Several distinct points are particularly 
convenient: 

• K*, the binding momenta of a trimer in the limit a ^ 00, 
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• a'^, the scattering length at which an Efimov trimer crosses the three-atom 
threshold, 

• a^, , the scattering length at which an Efimov trimer crosses the dimer-atom 
threshold . 

Other three-body parameters can of course be defined and generally all of them 
are related by simple expressions to each other, e.g. a* « 0.0798 [2]. 

A very different type of correlation plot can be generated by keeping the two- 
body scattering length constant while varying the three-body force. This allows 
the study of the correlation between different three-body observables. One such 
correlation that is very well known is the Phillips line (which will be shown later). 
It is an approximately linear correlation between the particle-dimer scattering 
length and the three-body binding energy. This correlation line is well-known in 
nuclear physics, since calculations of the neutron-deuteron scattering length and 
triton binding energy with different two-nucleon potentials lie near this line. 

Correlation lines calculated using the short-range EFT will in fact provide 
constraints on calculations with realistic potentials that give a large two-body 
scattering length. A three-body calculation of particle-dimer scattering length 
and binding energy employing such a potential has to give a point lying on or 
close to the Phillips line. The short-range EFT provides therefore strong low- 
energy constraints even in the absence of an experimental datum to fix the 
three-body low-energy constant. 

Higher Orders and Higher Partial Waves: 

The leading order (LO) calculations described above give results for the limit 
in which the range is taken to 0. The analysis of higher order corrections is 
important for several reasons. It allows to increase the accuracy of predictions 
for observables but provides also further information on the convergence radius 
of the short-range EFT. It was also mentioned before that the Lagrange density 
shown in Eq. ^ contains an infinite number of counterterms. It is therefore 
natural to ask at what order the next three-body force enters. 

The calculation of higher order corrections has been addressed a number of 
times. The correction to observables linear in the effective range was already 
considered by Efimov using the hyper spherical formalism [24^ [25] but no ex- 
plicit results were given for systems in which the Efimov effect is relevant. Ham- 
mer and Mehen pB] calculated the next-to-leading order (NLO) correction to 
neutron-deuteron scattering perturbatively. They found an improvement in the 
description of the corresponding phase shift but didn't present a result for the 
three-body binding energy at NLO. The extension of their perturbative approach 
to higher orders is involved since it requires the knowledge of the full off-shell 
three-body scattering amplitude. 

This and and other publications displayed an unexpected convergence pattern 
since the NLO shift to LO observables was unexpectedly small. This lead to a 
more detailed analysis of the NLO correction to observables which is linear in 
the effective range |27] . It was found that the exact discrete scale invariance of 
the leading order wave function in the unitary limit protects the LO bound state 
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spectrum. That means, if e'>P\^q (-^t^nlo) denotes the binding energy of the nth 
three-body bound state at LO (NLO) in the hmit a — > oo, then we have 

^tIo - ^t"nlo = ah n and 7 = . (27) 
This observation was made independently by Th0gersen, Fedorov and Jensen [28] 
using a numerical analysis of the Efimov bound state spectrum for finite range 
potentials. 

For arbitrary scattering length, we can define the linear shift in the binding 
energy of the nth three-body bound state in the following way: 



H/rp H/rp 



Fn[ — \+n,rsGn[^\ +0[{Krf] 



(28) 



and we have Gn(0) = for 7 = 0. In Ref. |27j the function G„ was also analyzed 
for finite 7 and evidence was found that Gn is close to at the atom-dimer 
threshold. An ongoing analysis will shine more light on the properties of the 
linear range correction and will address the impact of finite range corrections on 
universal relations that relate the recombination maxima at negative scattering 
length to the recombination minima at positive scattering length and will give a 
better understanding of how different signatures of Efimov physics are impacted 
by finite range effects [29|. 

Bedaque et al. [30j calculated the neutron-deuteron phase shifts up to next- 
to-next-to-leading order (N2L0) and included (based on a perturbative power 
counting argument that assumes the cutoff to be yl 1 /R) an additional energy- 
dependent three-body force at this order. They also suggested to include effective 
range corrections by modifying the two-body propagator in the STM equation. 
The two-body propagator with the effective range summed up to all orders is 
given by 

-^(^O'P) = , r , n / , ^ 77 • (29) 

-7 + -#(7 + n^Po - V IV + y-n^Po + P /4 



It is easy to show that the denominator in Eq. (j29p has two poles. In the case of 
positive 7 and this propagator has the bound state pole of the shallow dimer 
and a pole with binding energy 

16/1 1 



mrj \ 2 16 / 

This pole corresponds to a spurious bound state which is unphysical since in the 
two-body system the short-range EFT cannot make predictions for energy scales 
> l/{mR^). The spurious bound state pole in Eq. ([29]) leads to an ill-behaved 
kernel when inserted into the STM equation and the description of observables 
fails [31] . It was therefore suggested by Bedaque et al. to expand Eq. ([29|) in 
powers of the effective range: 

1 fs 'y^ + mpo — 



-D(po,p) 



-7 + x/-mpo + pV4 2 ^ ^_^pp + p2/4 
(-7 + V-mpo + P^/4 
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This amounts to a partial resummation of effective range corrections wlien in- 
serted into the STM equation. The difference in results for observables between a 
purely perturbative calculation of effective range corrections and this approach is 
expected to be of higher order. The solution of a modified STM integral equation 
leads directly to the desired observables and the amount of effort to calculate 
higher order corrections does therefore not increase compared to the solution of 
the leading order STM equation. 

In Ref. [32] a renormalization group analysis was performed for large cutoffs 
in the STM equation. It lead to the conclusion that for large cutoffs no additional 
three-body counterterm has to be included until N3L0 in the EFT expansion. 
Subsequent publications employing these results up N2L0 found good agreement 
with experimental measurements and calculations using realistic interactions [331 

There is an obvious disagreement between the predictions for the order at 
which the next three-body counterterm enters that requires a comment. In 
Ref. [30j it was argued that it enters at N2L0. In Ref. [H^ the RG analysis 
lead to the conclusion that N3L0 is the order at which an energy-dependent 
three-body force has to be included. The difference between the two analyses 
was the assumed size of the cutoff, namely either yl ~ l/r^ or yl ^ l/r^ which 
can in fact be the reason for the different conclusions. A future publication [29] 
will address this issue and will try to reconcile the results of Refs. [30] and [32] . 

An analysis of the power counting in higher partial waves was performed by 
Griei^hammer [35]. He defined the asymptotic exponent si which describes the 
momentum behavior of the half off-shell amplitude for angular momentum I and 
for large off-shell momenta k 

ti{k,p) oc fcV^' • (32) 

Using an argument based on the analysis of the perturbative evaluation of higher 
order corrections of the three-body amplitude, he identified the superficial degree 
of divergence of a given contribution at order n to be n — 2s;. A contribution at 
order n will therefore diverge if 

Re[n - 2si] > . (33) 

Naively, one expects the large momentum behaviour of this amplitude to be 
determined by the driving term Ki{k,p; E) of the STM equation 

ti{k,p) oc lim Ki{k,p,3k'^/A) oc k^/p^+^ , (34) 

which would imply the simplistic estimate si = I + 1. 

GrieBhammer derived, however, an exact algrebraic equation for the asymp- 
totic exponent si using a Mellin transformation of the integral equation defined 
in Eq. (12ip at zero energy 

l + s + 1 l-s + 1 2/ + 3 1" 
2 ' 2 ' 2 ' 4_ 

(35) 
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B^BG [mK] 
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127 


[2.186] 


125.5 


2.186 


4He4 


492 


128 


559.7 


132.7 



Table 1. Binding energies of the *He trimer and tetramer in inK. The two right columns show 
the Monte Carlo results by Blume and Greene [45] (denoted by the index BG) while the two left 
columns show the EFT results of Platter, Hammer and Meifiner [36]. The number in brackets 
was used as input to fix L3. 



It depends only on the relative angular momentum / and the function [a, b; c; x] 
is the hyper-geometric series. The order n at which the first three-body force in 
particular partial wave channel enters can then be obtained from a solution of 
Eq. (j35p and does for a number of cases not agree with the simplistic estimate 
derived from Eq. ([34]) . 

Equations (|35p and (|33p hold for spinless bosons. In Ref. [35] it was laid out 
how these formulas are generalized to systems of spin- 1/2 nucleons. 

2.4 The Four-Body System 

It is natural to ask whether a new counterterm has to be included for every new 
particle added to the problem. The question whether a four-body parameter is 
required for consistent renormalization in the four-body system is therefore a 
logical extension of the effort of applying the short-range EFT to few-body sys- 
tems. This pertinent question can also be paraphrased in more general terms: do 
zero-range two-body interactions and one three-body parameter lead to unique 
predictions for observables in the four-body sector? 

This question was addressed in the framework of the short-range EFT by 
Platter, Hammer and Meifiner in Ref. [36]. In this work the effective two- and 
three-body potentials at LO were generated and used together with the quantum 
mechanical few-body equations [37] to solve for the binding energies of the four- 
body system. The leading order effective two-body potential is given by 

(k|y|k') = X29{k)gik') , (36) 

here g{k) = exp(— fc^/yl^) is a regulator function. Since this is a separable poten- 
tial, the two-body problem can be solved exactly |38] 

tiE) = \g)TiE){g\ , (37) 

with the two-body propagator r 

T{E) = [l/X2-{g\Go{E)\g)r' , (38) 

and Go denotes the three-body propagator. Observables in the two-body sector 
will depend on the coupling strength A2 and the cutoff A. For a given cutoff 
A the coupling constant A2 is then renormalized by fixing the pole position of 




Figure 7. One-parameter correlations of three- and four-boson bound states. The crosses are re- 
sults for LM2M2 potential gS]. The triangles are the results for the TTY, HFD-B and HFDHE2 
potentials [?^H7] . 



the two-body propagator r which amounts to fixing the binding energy of the 
two-body bound state. The effective three-body potential is 

V3 = A3|0(ei, (39) 

where A3 denotes the three-body couphng constant that has to be adjusted to 
a three-body datum (such as the binding energy of the shallowest three-body 
bound state) and (U1U2IC) = exp(— (tif — 3u2/4)/yl^) is a regulator function (ui 
and U2 denote here the canonical Jacobi momenta in the three-body system |37j). 

Information on the necessity of a four-body parameter can be gained by 
studying the regulator dependence of four-body observables. A change in the 
regulator A corresponds to a modification of the short-distance physics while the 
renormalization conditions on the coupling constants A2 and A3 guarantee that 
low-energy observables in the two- and three-body sector remain unchanged. A 
detailed analysis shows then that no four-body parameter is required since the 
values of the tetrameter binding energies converge to well-defined values with 
increasing regulator A [36j. 

This approach was also used to address the four-body problem of "^He atoms. 
^He atoms have a scattering length 10 times larger than the range of the He-He 
interaction. The existence of n-body clusters of ^He atoms was shown experi- 
mentally [iQ], however, the measurement of the binding energies of these cluster 
states is currently not possible. Several potentials that are believed to describe 
the two-body interaction accurately have been developed and have been used in 
Faddeev [HI |42l [43l E] and Monte-Carlo [45] calculations. The four-body bind- 
ing energies were computed, using the results for the dimer and trimer binding 
binding energies obtained using the LM2M2 potential in Ref. [35] as input. For 
the three-body coupling constant, the excited three-body state of the ^He trimer 
was chosen as input parameter. The results can be seen in Table [TJ The bind- 
ing energies of the ^He tetramer were found to be in good agreement with the 
results for the ground and excited state of the ^He tetramer obtained by Blume 
and Greene [45]. 

By keeping the two-body parameter fixed and varying the three-body cou- 
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pling constant, the correlation line between the tetramer and trimer binding 
energies can be generated. These are are approximately linear correlations which 
are well known from calculations in the few-nucleon sector. They parameter- 
ize the common knowledge that an internucleon potential that gives the correct 
value for the triton binding energy also gives a good result for the four-nucleon 
bound state (i.e. the a-particle). In Fig. [7] we show this correlation line (which 
in nuclear physics is called the Tjon line) generated for values of the three-body 
binding energy close to the value of ^He trimer binding energy. The Tjon line is 
therefore a common feature of systems with a large two-body scattering length 
and does not depend on any details of the interaction at short-distances. 

This approach was used furthermore for a more detailed analysis of the four- 
boson system with large positive and large negative scattering length [39]. Results 
in this analysis also lead to the conclusion that every trimer state is tied to two 
universal tetramer states with binding energies related to the binding energy of 
the next shallower trimer: 

S4,o~5£'t and i?4^i ~ 1.01 E't for 7~0, (40) 

where E/j^^ denotes the binding energy of the deeper of the two tetramer states 
and -£^4^1 the shallower of the two. 

A recent calculation by von Stecher, d'Incao and Greene [H] supports the 
findings made in [36\ I39j . The authors of this work extended previous results 
to higher numerical accuracy. They furthermore considered the relation between 
universal three- and four-body bound states in the exact unitary limit (a — > cxd). 
They found 

^4,0 ~ 4.57^T and ^4,1 IMEt , (41) 

which agree with the results obtained in Ref. [39j and given in Eq. (j40p . 

The results obtained by the Hammer and Platter in Ref. [39] were furthermore 
presented in the form of an extended Efimov plot, shown in Fig. [8l Four-body 
states have to have a binding energy larger than the one of the deepest trimer 
state. The corresponding threshold is denoted by lower solid line in Fig. [8l The 
threshold for decay into the shallowest trimer state and an atom is indicated 
by the upper solid line. At positive scattering length, there are also scattering 
thresholds for scattering of two dimers and scattering of a dimer and two particles 
indicated by the dash-dotted and dashed lines, respectively. The vertical dotted 
line denotes infinite scattering length. A similar but extended version of this four- 
body Efimov plot was also presented by Stecher, d'Incao and Greene in Ref. [48]. 
They computed also the scattering lengths at which the binding energies of the 
tetramer states become zero and found 

a^Q « 0.43a* and a^i 0.92a* . (42) 

A further interesting detail was pointed out in [IS]. In the unitary limit the 
shallower of the two-four-body states can be considered as a trimer state with 
an extra particle attached to it. If this interpretation as an effective two-body 
state holds, then it leads to an important conclusion. Such an effective two-body 
picture in terms of a heavy and a light particle implies that the atom-trimer 




Figure 8. The a~^-K plane for the four-body problem. The circles and triangles indicate the 
four-body ground and excited state energies Bf^ and b'"^ \ while the lower (upper) solid lines 
give the thresholds for decay into a ground state (excited state) trimer and a particle. The 
dash-dotted (dashed) lines give the thresholds for decay into two dimers (a dimer and two 
particles). The vertical dotted line indicates infinite scattering length. All quantities are given 
in units of the three-body parameter L3. 



scattering length has to be large. For a two-body system with large positive 
scattering a length a and unequal masses, the binding energy is l/(2/xa^), where 
/X = 17111712/ {mi + 1172) is the reduced mass. With mi = m and m2 = 3m we 
obtain for cat 

Et - E,,o « , (43) 

3 ma^r^ 

where denotes the atom-trimer scattering length. We can therefore obtain a 
simple estimate for the scattering length of an Efimov trimer with index n and 
an atom in the unitary limit 



AT y 3m(i?T-i5^4,o) \/0.01m4") 

in) 

Since E}p is the only scale in the problem, naive dimensional analysis predicts 
the range of the atom-trimer interaction to be of order The estimated scat- 
tering length OAT is thus by an order of magnitude larger than this estimate of 
the range. 



3 Lovif-Energy Universality in Atomic Physics 

In this section we will discuss recent applications of the short-range EFT to 
systems of cold atoms. We will focus in particular on the problems of three-body 
recombination and atom-dimer relaxation in ultracold gases but will also report 
on recent progress in the four-body sector. 
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The recombination rate and atom-dimer relaxation have been identified as 
key signatures in the search for Efimov physics since they can be obtained my 
measuring atomic loss rates. Three-body recombination is a process in which a 
two-body bound state is formed as the result of a three-body collision. The two 
outgoing particles will gain kinetic energy in this process equal to the binding 
energy of the dimer. In experiments with ultracold atoms, the kinetic energy is 
often sufficient to allow the atom and dimer to subsequently escape the trapping 
potential. 

Ultracold gases of alkali atoms are very well suited for such experiments since 
in many cases Feshbach resonances allow the scattering length to be tuned to ar- 
bitrarily large values using an external magnetic field. Observables can therefore 
be measured as a function of the two-body scattering length and results can be 
compared to the corresponding one-parameter correlation discussed previously. 

The first experimental evidence for Efimov physics in an atomic system was 
presented by Grimm and co-workers [l9]. In this pioneering experiment with 
ultracold ^^^Cs atoms in the lowest hyperfine state, they observed a resonant 
enhancement in the 3-body recombination rate at a ~ —850 ao that can be 
explained with an Efimov trimer close to the 3-atom threshold. 

Since then the number of experiments that show evidence for universal three- 
body physics has increased significantly. Signatures of Efimov physics have been 
found in three-component Fermi gases of ^Li \50\ I51j . in a Bose gase of 
atoms [52], and in heterogenuous mixtures of "^^K and ^^Rb [53] • The level of 
sophistication in experiments has increased so much that recent measurements 
even test the implications of low-energy universality on four-body dynamics |54j . 

One experimental tool that has become essential for the analysis of Efimov 
physics are Feshbach resonances. As mentioned above they allow the interparticle 
scattering length to be tuned by adjusting an external magnetic field. A Fesh- 
bach resonance arises due to the coupling of two atoms in an open channel to a 
closed channel. The open channel corresponds to a pair of atoms in energetically 
allowed hyperfine states while the closed channel corresponds to a combination of 
hyperfine states that is energetically inaccessible to asymptotic scattering states. 
The coupling between the open and closed channel arises due to the hyperfine 
interaction. The scattering length becomes large when a magnetic field is used 
to tune a bound state in the closed channel to the threshold of the open channel. 
The dependence of the scattering length on the magnetic field near a Feshbach 
resonance can be described by 



where Bq denotes the scattering at which a state in the closed channel is at 
threshold and A governs the width of the resonance. For magnetic fields away 
from the resonance Bq, the scattering length is given by Obg. The scattering 
length can therefore be tuned to values significantly larger than the low-energy 
length scale of the atom-atom interaction which is given by the so-called van der 
Waals length iydW- This quantity is related to the long-range van der Waals tail 




(45) 
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of the atom-atom interaction via 

V{r) ^ -^^^^ . (46) 

This length scale sets therefore the natural scale for the effective range Vg ^ R 
of the interaction and quantifies at which distances the short-range EFT is not 
applicable anymore. 

3.1 Thrcc-Body Recombination of Identical Bosons 

In a three-body recombination process, three particles collide and two of them 
form a two-body bound state. If the scattering length is large and positive, the 
resulting two-body bound state can be a deep dimer or a shallow dimer when the 
scattering length is positive. A deep dimer is a two-body bound state that cannot 
be described within the short-range EFT and has binding energy > l/(mR?). 
The shallow dimer can be described within the short-range EFT and has binding 
energy ~ l/(ma^). If the scattering length is large and negative the resulting 
two-body bound state can only be a deep dimer. 

The change in density n due to such losses is described by 

A„ = -niost" n , (47) 

where niost« = -^3 is an experimentally measurable loss rate constant, niost is the 
number of atoms escaping in every recombination process and a is the event rate 
of the recombination process. In the case of a Boltzmann distribution, the event 
rate can be related to the hyperangular averaged recombination rate K(E): 

"^^^ ~ 6 J^-' dE E^ e-mf^sT) ■ (48) 

The recombination rate K{E) can be decomposed into a contribution from re- 
combination into the shallow dimer and a contribution from recombination into 
deep dimers: 

K{E) = i^shallow(^) + ^deep(^) • (49) 

The three-body recombination rate into the shallow dimer can be decomposed 
into channel contributions with different total orbital angular momentum J: 

oo 

K,^sIIoAE) = Y,k(J\E) . (50) 

J=0 

The recombination rate -f^shallow into the shallow dimer is then related to the 
total atom-dimer breakup cross section via 

-K'shallow(-E) = ^ J^2^2 ^ ^ '^breakup(-E^) • (51) 

This implies that the recombination rate is related to the S-matrix for the scat- 
tering of three atoms into an atom and dimer through the following relation: 

^<^'(^)-^^^^Si^Ei4ti,..(^)r. (52) 

n=3 
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Here, n denotes a set of quantum numbers that includes the relative angular 
momenta between the particles. The hyperangular average is implemented by 
the sum over n which starts at n = 3 for convenience. The unitarity of the 
S-matrix in the angular momentum J sector implies 



s'xiADm' + E &AAAiE)\' = 1- (53) 

n=3 



By unitary, the recombination rate into the shallow dimer is therefore directly 
related to the S-matrix element for elastic atom-dimer scattering: 



(54) 



smce 



S^AADiE) = e'''--^'^\ (55) 



The phase shifts for atom-dimer scattering are therefore sufficient for a calcula- 
tion of the recombination rate into shallow dimers. 

Recombination into the Shallow Dimer: 

Bedaque, Braaten and Hammer were the first to use the short-range EFT to 
describe three-body recombination into the shallow dimer [55]. They calculated 
the three-body recombination rate at zero temperature and for J = by solving 
the integral equation for the scattering of three atoms into an atom and dimer, 
shown diagrammatically in Fig. [9l It can be shown that for zero energy, the 
amplitude for this process is related to the amplitude for elastic scattering given 
in Eq. (j22|) at an off-shell point. This can be used to simplify the numerical 
calculation of the three-body recombination rate. This approach can also be 
used to calculate the contribution to the recombination rate from channels with 
J 7^ 0. The contribution to the three-body recombination rate is in this case 

J,,r,^E,^iJ±^3f^±MlM'l, (56) 

m 

where fj{x) is a real- valued scaling function: 

/j(x) = l-exp(-4Im5i5,(^)), (57) 

and X is defined as 

x={ma^Ef''^. (58) 

As X ^ 0, the leading powers of x are determined by Wigner's threshold law 
[56]: fj{x) ~ where Ai = 3 and \j = J for J > 2. 
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It is only in the channel with total angular momentum J = that observables 
depend on the three-body parameter. In this channel, a different perspective on 
the recombination rate can be gained using Efimov's radial [641 [2] laws which 
expresses the S-matrix in the three-body sector as a combination of universal 
functions of a scaling variable x and a phase 6^q: 

(j=o) .p. . S2i(aOV*^ 

^ADADiE) = ^22(x) + ^— ^^^^^y^ , (59a) 



(J=0,n) .p. ^ ^ , ■S2l(x)gln(x)e^^^*° 

SAD,AAAiE) = S2n{x)+ ^ _ (^)g2ie,o " ^^^^^ 



The phase 0,,o is related to the minimum in the recombination rate a*o for positive 
scattering length 

6**0 = soln(a/a*o) , (60) 

where a*o ~ [2]. 

Using analytical results on the atom-dimer phase shift at breakup threshold 
obtained in Ref. |57j, Efimov's radial laws can be used to derive an analytical 
result for the recombination rate into the shallow dimer 

j^(o)(^ = 0) = ''^^^^(^^ ~ 3\/3)sin^[soln(a/a^o)] a'^ ^^^^ 
sinh^(7rso) -|- cos^[so ln(a/a*o)] 

Macek, Ovchnikov and Gasaneo derived this analytic result for the three-body 
recombination rate at the three-atom threshold in [58]. This equation was also 
derived independently by Petrov [59j. Simpler approximate expressions that cor- 
respond essentially to omitting the cos^ term in the denominator have been 
previously derived by Nielsen and Macek [60], Esry, Green and Burke [61] and 
by Bedaque, Braaten and Hammer [55] . 

For a fixed scattering length, the recombination rate in Eq. (j6ip has therefore 
the maximum value 

-K'max = 6Cmaxa^M with Cmax = ^^^^^-^^7 (62) 

smh [ttsq) 



Effects of Deep Dimers: 

While recombination into the shallow dimer can only occur for positive scatter- 
ing length, recombination into deeply bound states (into deep dimers) can occur 
for either sign of the scattering length. Properties of these bound states can- 
not be calculated in the short-range EFT. However, their cumulative effects on 
low-energy observables can be accounted for analytically continuing the Efimov 
parameter to complex values. This introduces only one additional real- valued 
parameter denoted usually with ry^,. If the dependence of an amplitude on or, 
equivalently, a*o is known analytically, the effect of deep dimers can be taken 
into account by the simple substitution 



lna*o — > lna*o — irj*/so 



(63) 
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For positive scattering length, making this replacement in the amplitudes for 
recombination into the shallow dimer leads to 

7687r2(47r-3^/3)(sin2[soln(a/a*o)] +sinh2r?,) 

smh (vrso + + cos^[so ln(a/a^,ojJ 

Using unitarity, one can also obtain an analytic expression for the recombination 
rate into deep dimers: 

3847r^(47r - 3^3) coth(^3o) sinh(27?,) 

^deep(Uj = ——21 \ T- ^7 -r~!~l ' l« > Uj • (05 j 

smh (vrso + ?7*) + cos^[soln(a/a^,o)J 

These results were first derived by Braaten and Hammer [2|. Simpler approx- 
imate expressions that correspond essentially to omitting the cos^ in the de- 
nominator were derived in Ref. [62] • The weak dependence on the three-body 
parameter a*o in Eq- (|65|) was first observed in a numerical calculation for the 
case of infinitesimal ry* [63] . 

For negative scattering length , there is no shallow dimer. The recombination 
rate into deep dimers is given by 

/n^ 765sinh(27?^) 

-?^deep(0) = .2/ w / /^^ I -2 ' (« < 0) • (66) 

sm (so ln(a/a'^)) + sm m 

This results was first derived by Braaten and Hammer in Ref. The scaling 
of -ftTdeep with was predicted by Nielsen and Maceck and by Esry et al. [60l [61] . 
Esry et al. were the first to point out the existence of a log-periodic sequence of 
resonances related to Efimov trimers [61] • 

Three-Body Recombination at Finite Temperature: 

The experimental measurements of the three-body recombination rate of ^^^Cs 
atoms were performed at finite temperature. Although the zero temperature 
equations can be applied immediately to their results, it is desirable to extend 
these to finite temperature. This requires the calculation of the energy depen- 
dent recombination rate. For recombination into the shallow dimer, this can be 
obtained by solving the amplitude for scattering of three atoms into atom and 
dimer as a function of the external momenta, calculating the thermal average of 
the recombination rate, and then adjusting the three-body parameters and r/^, 
to fit experimental data. An alternative path is to determine the universal func- 
tions Sij{x) from the atom-dimer scattering phase shifts. Once these functions 
are determined, the temperature dependent recombination rate can be deter- 
mined for any system of identical bosons with arbitrary three-body parameter. 
This was done in Ref. [65] for the case of positive scattering length. Using the 
STM equation the atom-dimer phase shifts were calculated for a wide range of 
energies and three-body parameters. The results were then fit to the formula 

/ . (j=o).ps\ , ^ , sujxf exp[2iso ln(a/a^o)] 

exp[2i5\j^ '{E}] = S22ix) + - 1—^ t^. — t-^, r, 67 

V / 1 - sii(x) exp[22So ln(a/a*o)J 

to determine the universal scaling functions sii(x), suix) and S22{x) for < 
X < 10. 
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Using the unitarity of the S-matrix it can then be shown that these three 
universal scahng functions are sufficient for the calculation of the J = con- 
tributions to the recombination rate into shallow dimers, even if there are deep 
dimers: 



144\/37r2 



1 



2 2ie,o-2v* 



1 - Sii(x)e2»e*0-2r?, 



Sl2(X 



a 
m 



|1 - sii(x)e2*^*o-2^*|2 
They are also sufficient to calculate the recombination rate into deep dimers: 



(68) 
(69) 



144^/37r^(l - e-^^'){l - |gii(x)p - ^x^) 



(70) 



The knowledge of three universal scaling functions facilitates therefore the calcu- 
lation of the temperature dependent recombination rate constant in the presence 
of deep dimers. In Fig. [10] the theoretical results for three different values of r/* 
are compared to the experimental results by Grimm and co-workers. The vertical 
axis is the recombination length defined by 



P3 



2m ^ 



(71) 



where niost- 

The recombination minimum that determines the three-body body input is 
near 200 oq. This is comparable in size to the Waals length scale (mCe)^/^ ~ 
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200 oo, so range corrections may be large near the minimum and the disagreement 
between theoretical and experimental results at smaller values of the scattering 
length is hardly surprising. The authors of Ref. [65] were not able to determine 
r/^, since their fit was insensitive to the value of ij^,. It yielded, however, the 
upper bound rj^ < 0.2. Their result for the recombination minimum amin = 
210(10) ao obtained from a fit for a > 500 oq agrees with direct experimental 
loss measurements preformed by the Innsbruck group. 

Finite Range Effects in Three-Body Recombination: 

The recombination rate measurements performed by Grimm and co-workers 
were carried out at scattering lengths at which the finite range of the atom- 
atom interaction is expected to play a significant role. The range of the ^^'^Cs 
interaction is of the order of 200 ao- It is therefore of interest to understand 
the impact of the finite range of the atom-atom interaction on experimentally 
measurable quantities. 

A first effort to calculate range corrections to the recombination rate into 
shallow dimers was performed in [66]. The authors considered corrections up 
to N2L0 in the EFT expansion to the recombination rate and calculated the 
correlation between the atom-dimer scattering length and the recombination 
rate K{0) at zero energy. The approach was furthermore used to calculate the 
recombination rate coefficient for ^He atoms for which the effective range is 
known. This approach can be applied to more interesting systems such as ^^^Cs 
or ^Li provided the corresponding effective range is known. 

3.2 Three-Body Recombination of Fermions 

A large number of experiments with cold atoms are now carried out with fermionic 
atoms. Systems of fermionic atoms are of interest due to their relation to solid 
state physics. In particular, systems with two spin states have received con- 
siderable attention. Many-body systems of cold atoms display superfiuidity at 
sufficiently low temperatures. Feshbach resonances can be used to study how the 
mechanism for superfiuidity depends on the interactions. In the case of fermions 
with two spin states, as the scattering length is varied from l/a<Otol/a>0 
through the Feshbach resonance, the mechanism changes continuously from the 
formation of Cooper pairs to the Bose-Einstein condensation of shallow dimers. 

The Efimov can not occur in systems of fermions with only two spin states. 
This can be understood from the fact that a pointlike three-body S'-wave interac- 
tion is forbidden due to the Pauli principle. However, in the case of fermions with 
three spin states the Efimov effect can occur again since a pointlike three-body 
5-wave interaction is not forbidden by the Pauli principle. Experiments using ^Li 
atoms in the three lowest hyperfine states have been performed by Selim and co- 
workers [50] and O'Hara and co-workers [51]. The three corresponding scattering 
lengths in this system have broad Feshbach resonances at nearby values of the 
magnetic field, i.e. at 690 G, 810 G and 833 G. A narrow loss feature near 130 G 
was discovered independently by both groups. At these magnetic field strengths 
all three scattering lengths are negative and relatively large, universal results are 
therefore relevant in this region and the short-range EFT is applicable. 
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Figure 11. The 3-body recombination rate constant Kz for the three lowest hyperfine spin 
states of ®Li atoms as a function of the magnetic field B. The two vertical dotted lines mark 
the boundaries of the region in which |oi2| > 2 l^dw- The solid squares are data points from 
Ref. [50]. The solid dots are data points from Ref. [51]. The curve is a 2-parameter fit to the 
shape of the data from Ref. [50] . 



The rate equations for the number densities of atoms in the three spin 
states are 

d 

—Ui = -K^nin2n^. (72) 

In the case of this more comphcated system, no analytical results for the re- 
combination rate are available so one has to do a numerical calculation. The 
three-body recombination is related to the amplitudes for elastic scattering 

= V a,ajWy (0,0), (73) 

The amplitudes Aij have to be calculated from a set of nine coupled integral 
equations which are an extension of the STM equation: 



-4(g,0) 

A ~ . 



(74) 



. -4(g,0) 
■ -l/ofc + ^/3/2g ' 

These equations are written out without the use of an explicit three-body force 
and instead the momentum cutoff is used to fix the the three-body coefficient. 
In |67j the recombination rate of general systems of three identical fermions with 
three spin states was analyzed. The case of ^Li was also considered. Theoretical 
results on the magnetic field dependence of the three scattering lengths were 
used to calculate the three-body recombination rate in regions of the magnetic 
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field relevant to the Penn-State and and Heidelberg experiments. The effects of 
deep dimers were included by using a complex cutoff Ae'^'^*^^^ in the STM integral 
equation in Eq. ([71]) . The three-body parameters were adjusted to reproduce the 
observed recombination maximum near 130 G. In Fig. [11] we show a comparison 
of theoretical and experimental results. The results are in excellent agreement 
with the observed recombination rate near the narrow loss feature at 210 G. 
They also predict a second narrow loss feature near 500 G. In the experiments 
a broad loss feature is observed near 500 G, but the behavior in this region is 
not correctly reproduced by theory. The reason for this discrepancy is under 
investigation. 

Following this work, the problem was also considered using a wave function 
approach [68] and functional renormalization [69j. The authors of both studies 
found qualitative agreement with the EFT results. 



3.3 Atom-Dimer Relaxation 

Another observable closely related to three-body recombination is atom-dimer 
relaxation. If there is an Efimov trimer close to the three-atom threshold, inelas- 
tic scattering processes of atom-dimer to atom- (deep dimer) will be resonantly 
enhanced. The parameters that determine the rate of these losses are a, and 

The relaxation rate event constant is defined by the following equation: 

= ^riD = -PnATiD , (75) 
at at 

where ua denotes the atom density and denotes the density of dimers. The 
rate constant /? can also be obtained using the optical theorem: 

P = -— Im OAD , (76) 
m 

where oad denotes the atom-dimer scattering length: 

20.3sinh(2r?,) a ^^^^ 

sin^(so ln(a/a*)) + sinh^ r/* m 

where a^: denotes the value of the two-body scattering length at which there is a 
trimer at the atom-dimer threshold: a* ~ 0.0798 k~^. This expression was first 
derived in |62j using Efimov's radial law discussed in the previous section. 

Atom-dimer relaxation at finite temperature was first considered by Braaten 
and Hammer in Ref. [70]. Recently, Helfrich and Hammer [TTj extended this 
work. Using generalized Bose-Einstein distribution functions to perform a ther- 
mal average, they calculated the rate of change in the number of shallow dimers 
Nd due to relaxation into deep dimers: 

d /" 3k 

-Nd = - J d^rd^pAd^PDnAiPA,r)nDiPD,r)—aTD > (78) 
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Figure 12. Left: The dimer relaxation coefficient /3 as a function of a/ao for T = 170 uK, 
a, = 395ao, and different values of »7, . BO indicates a Boltzmann average. Right: /? as a 
function of a/ao for T = 40 nK, a* = 395ao, and different values of r;*. The data points in both 
panels are from [72] . 



here k = |2p^ — Pd|/3 denotes the relative momentum, p^i and po are the 
momenta of atoms and dimer, respectively. The generalized Bose-Einstein dis- 
tribution functions ny^/^j for atoms/dimers in a harmonic trapping potential are 



niipi,r) 



exp 



Pi 



+ 



2 2 

mjijo r 



2mi 2 



1 



(79) 



with i = A,D, oj denotes the average trap frequency, and /Xj is the chemical 
potential. Knowledge of the 5-wave phase shift for atom-dimer scattering 



allows therefore the evaluation of Eq. ([78]) since a 
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(80a) 



Braaten and Hammer have calculated the phase shift ka cot 6q^ for k up to the 
dimer breakup threshold \/3/(2a) [62j. They parameterized the results as 



kacotdo^ik) = ci{ka) + C2ika) cot[soln{0.19 a/a^) + (p{ka)] , (81) 



where 



ci(fea) 
C2{ka) 
4>{ka) 



-0.22 + O.SO/c^a^ 
0.32 + 0.82fc2a2 _ 



O.lTfc^a^ 



2.64 



QMk'^a^ + 0.23^^0"^ , 



(82) 



Using these results Helfrich and Hammer evaluated the integral in Eq. (I78p . They 
related there results to the experimental results by Knoop et al. [72] using the 
loss model employed in [72]. They find 
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(83) 
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The factor of (\/3/2)'^ takes into account the assumed Boltzmann distributions 
exp(— mja;^r^/3) of the atoms and the dimers that were assumed in the analysis 
of Ref. [72]. 

The free parameters a* and r/* are fitted to the data at T = 170 nK. In the 
left panel of Fig. [121 the results for a* = 395ao and rj^ = 0.025, 0.06 (solid curves) 
give an excellent fit to the experimental results by Knoop et al. [72\. The dashed 
curve in the left panel is the result one obtains by using a Boltzmann instead of 
a Bose-Einstein distribution in Eq. (j79p . The results display some sensitivity to 
r]^. At T = 40 nK, the best fit is obtained with 7]^, = 0.012 (solid line). The value 
r]^ = 0.025, which gave the best fit for 170 nK gives a prediction that is too large 
by a factor of two when compared to the experimental results. 

3.4 Four-Body Recombination 

Stecher, d'Incao and Greene considered in their work on the four-body system 
also the implications of 4-body universality on loss rates in systems of ultracold 
atoms [73j. In particular, they pointed out that the existence of the two universal 
tetramer states discovered by Platter and Hammer should lead to observables 
loss features at values of the scattering length that are related to the scattering 
length at which the tetramer binding energy becomes zero. Their results for the 
values of the scattering lengths are given in Eq. (|42p . 

In a second paper Stecher, d'Incao and Greene addressed the problem of 
dimer-dimer collisions [73] which is relevant for the case of positive scattering 
length. They denoted the two-body scattering length for which the ith four- 
body bound state crosses the dimer-dimer threshold by a^di- They denoted the 
two-body scattering length at which the trimer-atom and dimer-dimer channels 
become degenerate by a'^. We previously defined o^, to be the scattering length at 
which the trimer state crosses atom-dimer threshold. The ratio of these quantities 
are universal numbers 

«dd,o/«dd ~ 2.37 , aJd i/a;^^ ~ 6.6 , and a^^Ja^ ^ 6.73 . (84) 

Ferlaino et al. recently studied the four-body problem with short-range interac- 
tions experimentally [51]. Using ultracold ^^^Cs atoms in the lowest hyperfine 
state at a temperature of 50 nK, they found loss features at scattering lengths 
— 730ao and — 410ao which were interpreted as the four-body loss features pre- 
dicted by Stecher, d'Incao and Greene [l8]. With the triatomic Efimov resonance 
measured at — 870ao, this gives for the ratios of the four- and three-body reso- 
nance position 

04 o/a* PS 0.47 and 04 ^/a* w 0.84 . (85) 

These experimental results are in fact surprisingly close to the zero-range pre- 
diction made in [18] since finite range effects are expected to important at these 
values of the scattering length. The range of the Cs-Cs interaction (which is set 
by the van-der Waals length scale) is approximately 200ao. 
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3.5 Challenges and Opportunities I 

Recent experimental progress demonstrates that the hmits of complexity have 
not been reached yet in the field of few-body dynamics in gases of ultracold 
atoms. The experimental evidence for three-body universality in systems of ultra- 
cold atoms has increased significantly. First experimental evidence for four-body 
universality have been presented and more can expected in the near future. 

Measurements of Efimov loss features in heterogenous systems, i.e. systems 
with different constituent masses, have also been presented. It is therefore a 
natural task to identify the universal signatures and their relation to each other 
using the short-range EFT. Some of these signatures might offer also stronger 
signatures of Efimov physics since the discrete scaling factor associated with the 
Efimov bound state spectrum depends on the mass ratios of the constituents. 

Significant progress in mapping out the relevant signatures of universal four- 
body physics has been made. However, a full treatment of scattering processes 
remains desirable. This would facilitate an analysis of the impact of universality 
on further quantities such as the dimer-dimer or atom-trimer scattering lengths. 
Access to scattering quantities is also required for the calculation of the four-body 
recombination rate at finite temperature. 

The ongoing analysis of the effects of range corrections on universal few- 
body physics is important since most experiments are carried out at scattering 
lengths at which the finite range of the two-body interaction is expected to have 
a measurable effect. The questions that remain to be answered are therefore how 
universal relations that relate different recombination observables will be affected 
by finite range effects and what implications does a small but finite range have 
for the four-body bound state spectrum. Since the linear range correction has 
been shown to be small for some observables it will be important to address the 
quadratic range correction with a full N2L0 analysis. 

4 Low-Energy Universality in Nuclear Physics 

The model-independent description of nuclear systems has been one of the main 
goals in nuclear physics for many years. Since their introduction to nuclear 
physics, EFTs have allowed the calculation of a large number of observables 
model-independently and in some cases to very high accuracy. The standard 
EFT approaches to internucleon interactions employ nucleon and pion fields as 
the minimal set of degrees of freedom |741 HI [5] . The separation of scales that 
is exploited here is the one of chiral symmetry which identifies pions as the 
Goldstone bosons of the spontaneously broken chiral symmetry of QCD. The 
relatively small pion mass is then a consequence of the small up and down quark 
mass, i.e. the explicitly broken chiral symmetry. 

Fortunately, this is not the only separation of scales in the two-nucleon sys- 
tem. It turns out that also the scattering length in the nucleon-nucleon system 
is large compared to the range of the interaction which in this case is set by 
the pion mass m^. The short-range EFT can therefore be applied in the nuclear 
sector as long as we consider momenta that are much smaller than ttZtt. 

Such a limitation might seem to make this approach inadequate for nuclear 
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physics. However, it is important to realize that a large number of reactions 
relevant to nuclear astrophysics occur at energies well below this breakdown scale 
of the short-range EFT. Proton-proton fusion {p+p — > ^H+e"*" + fe) for example, 
a reaction which plays an important role in the sun's energy generation, occurs 
at the keV scale due to the low temperature in the sun. The short-range EFT 
seems therefore to be the ideal framework for the calculation of such reaction. 

Halo nuclei and weakly bound systems of a particles are an additional play- 
ground for the short-range EFT. Halo nuclei are weakly bound systems of a core 
(e.g. an a-particle) and additional nucleons that are weakly bound to the core. 
There is also evidence that some of the properties of a-clusters can be described 
using the short-range EFT. 

In this section we will discuss applications of the short-range EFT to systems 
of nucleons. We will first consider the two-nucleon system and discuss a recent 
calculation in this sector which exemplifies the strong predictive power of the 
short-range EFT. We will then turn to the three-nucleon sector. We will discuss 
the renormalization, the inclusion of higher order corrections and the status of 
the inclusion of external currents in the few-body sector. Then we summarize 
recent calculation of halo nuclei. We end with a subsection which discusses the 
work required to be done in the future. 



4.1 Two Nucleons 

Nucleons are spin 1/2 particles that come in two flavors, protons and neutrons. A 
nucleon field will therefore carry spin and isospin indices to accommodate these 
extra degrees of freedom. The Lagrange density for nucleons interacting at very 
low energies is constructed by writing down all possible operators allowed by the 
underlying symmetries 

-C^ {N^a2TaT2Ny {N^a2TaT2N) + ... , (86) 

where the dots represent higher-order contributions suppressed by more fields 
and/or derivatives. The low-energy constants Cg and Cq are renormalized to the 
spin-singlet and triplet scattering length, respectively. 

The short-range EFT has been applied very successfully in the two-nucleon 
sector to electroweak observables. Here, we will only mention a few very recent 
calculations since they exemplify typical applications and the predictive power 
for low-energy processes with only a few physical input parameters. 

Christlmeier and GrieBhammer [75] recently addressed a discrepancy be- 
tween previously obtained theoretical and experimental results for the electro- 
disintegration cross section of deuterium (^H(e, e'p)n). They considered the triple- 
differential cross-section for this process which can also be decomposed into 



+aLTCOs0p + aTTCOs2^p) , (87) 
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where i^l^^ = {0^^^,<l'^^^ = 0) and Qp = {0p,^p) are the scattering angles of 
electron and proton in the lab frame, respectively. The superscript "lab" denotes 
the laboratory frame which is defined to be the rest frame of the deuteron. 
They found their results to be in excellent agreement with previous theoretical 
calculations by Arenhovel et al. [76] and Tamae [77] which lead to the conclusion 
that none of the available theoretical approaches can explain the experimental 
data presented for glt in Ref. f78]. A subsequent experimental study of the 
double differential cross section for (^H(e, e'p)n) at an angle of 180° lead to 
excellent agreement between experiment and theory [79j . 

Finally, Ando considered pp ppir^ near production threshold in the short- 
range EFT [80j and Phillips, Schindler and Springer considered parity violating 
nucleon-nucleon scattering in |81j . 

4.2 Three Nucleons 

The short-range EFT was first applied to three-nucleon systems by Bedaque, 
Hammer and van Kolck [82^ I83j. As in the bosonic case, the calculation of three- 
body processes is simplified if we introduce auxiliary fields that carry the quan- 
tum numbers of the allowed two-body 5-wave states. Here, we introduce two 
auxiliary fields (one for every possible S-wave channel), ti with spin 1 (isospin 
0) and Si with spin (isospin 1), respectively 

^ = N^iido + ^)N-tl{ido + ^-At)U-s}{ido + ^-As)s, 



+gt [tlN' T2<Ti(T2N + h.cj + Qs (^sjiV' <T2<TjT2N + h.C 

-G^N^ gj{tiai)himii + \gtgs[{ti(yi)^ SjTj + h.c] 



+gl{SjTj)hjiTj^N + . . . 



We want to use this Lagrange density to calculate observables for nucleon 
deuteron scattering. The total spin of this system is either 1/2 (singlet) or 3/2 
(quartet) and correspondingly two different integral equations can be derived 
that describe the scattering of a neutron and a deuteron in a relative S-wave. 
The integral equation for scattering in the quartet channel is given by 

h/2{p,k) = -K{p,k) - - / dqq''Ds{q)K{p,q)t3/2{q,k) , (89) 

where 

1 ( p^ +pk + k'^ 



The Pauli principle forbids the appearance of an S-wave three-body force in this 
channel, since the spins and isospins of two of the three nucleons are aligned 
parallely. 
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The integral for the scattering in the singlet channel is given by 

+ dq^Dt{q){K{q,p) + ^%{q,k) 

A/2{P^k) = ^(^^K{p,k) + ^^+^j\q'Dt{q){K{p,q) + ^)ts{q,k) 

rA n TT 

■ dq^Ds{q)iKip,q) + ^)tsiq,k) . 



(91) 



In this case a three-body force is allowed and also required for consistent renor- 
malization of the problem. The triton binding energy can be calculated from the 
homogeneous part of Eq. (I9ip since the total spin (isospin) of the triton is 1/2 
(1/2). 

The solutions of Eqs. ([89]) and (|9T|l are related to the deuteron-neutron phase 
shift 

Svr 1 

t\,Jk,k) = —- -. (92b) 

' ^ m A;cot(5i/2 - ik ^ ' 

The appearance of a three-body force in the doublet channel implies the corre- 
lation of different three-body observables. The short-range EFT describes these 
correlations therefore with the minimal number of degrees of freedom and all ob- 
served correlation lines are therefore a consequence of the large scattering length 
in the two-nucleon system. The short-range EFT also offers a different perspec- 
tive on the internucleon interaction: different two-nucleon interactions may de- 
scribe two-body data equally well but give different predictions for three-body 
observables. 



Finite Range Corrections: 

Since the ratio of effective range over scattering length in the spin-triplet 
channel is approximately R/a ~ 1/3, the consistent calculation of finite range 
corrections is crucial for an accurate description of low-energy observables in nu- 
clear physics. In the three-nucleon context, range corrections were first addressed 
by Hammer and Mehen [26]. They included the the first subleading correction 
perturbatively. Bedaque et al. [30] calculated finite range corrections up N2L0 
in the EFT expansion using the modified version of the STM equation discussed 
in subsection 12.31 They calculated phase shifts for neutron-deuteron scattering 
up to N2L0 using the triton binding energy and the neutron-deuteron scattering 
length as three-body input parameters. In [33] the powercounting results from 
Ref. [32] were used to calculate scattering phase shifts for 5-wave scattering in 
the neutron-deuteron doublet channel up to N2L0. The triton binding energy 
was calculated as a function of the neutron-deuteron scattering length 03. In 



36 



Low-Energy Universality in Atomic and Nuclear Physics 




Fig. [T3]we show the result for this correlation line (the Phillips line). The dot- 
dashed, dashed and solid line show the LO, NLO and N2L0 results, respectively. 
The cross denotes the experiment result for the triton binding energy and the 
neutron-deuteron scattering length. The lines in Fig. [13] are not parallel to each 
other since the shift from order to order depends on two expansion parameters, 
KVg and (with k = ^JmBt). At the experimental value of the nd scattering 
length we find for the triton binding energy 

Bt = (8.08(LO) + O.ll(NLO) + 0.35(N2LO)) MeV . (93) 

These 8.54 MeV are in very good agreement with the experimental value Bf^^ = 
8.48 MeV and within the naively expected error bounds of a N2L0 calcula- 
tions. These results furthermore exemplify the unusual convergence pattern of 
the short-range EFT for three-body observables. This is due to the unnaturally 
small NLO correction while the N2L0 corrections has the size as predicted by 
naive dimensional analysis. We explained in Section [2.31 that in the unitary limit 
the discrete scale invariance of the LO wave function prohibits the bound state 
spectrum to obtain any NLO correction in the unitary limit. It is tempting to 
speculate that the approximate discrete scale invariance of the leading order am- 
plitude at unitarity is also partially responsible for the unexpected convergence 
pattern in the three- nucleon sector. 

Electroweak Observables: 

A major benefit of using a field-theoretic framework is the straightforward 
inclusion of external currents. The short-range EFT facilitates therefore the cal- 
culation of electroweak observables of few- nucleon reactions. This is particularly 
useful since a number of reactions relevant to nuclear astrophysics occur at en- 
ergies well below the breakdown scale of the short-range EFT which is set by 



L. Platter 37 




7 8 9 

[MeV] 

Figure 14. The correlation between the triton charge radius and binding energy. The solid 
(dashed) line denotes the leading-order result using at and Oa {Bd and as) as input parameters. 
The circles indicate Faddeev calculations using different potentials from Ref. [87] while the 
square gives the experimental values. 



the pion mass. The number of three-body calculations with external currents is 
however, extremely limited. Universal properties of the electric form factor and 
charge radius of the triton were considered in [8l] . The leading order wave func- 
tion was used to evaluate the charge form factor from the leading order charge 
operator 

Fc{^^) = {^K+d k/|pc|ki i^^k) , (94) 

where q = kj — kj, kj and kj are the initial and final momentum of the scat- 
tered electron, and I^'k denotes the full triton wave function with center of mass 
momentum K. The charge density operator pc is defined as ^85j 

3 . 
PC = T.y^+n.)p'cir-r,) + -il-n,)pUr-r,)] . (95) 

i 

The charge radius (r^) can then be defined as 

Fc(q^) = l-q^(r2)/6 + ... . (96) 

In Fig. [13] we display the one-parameter correlation between the charge radius 
and the triton binding energy. The two lines shown in this figure correspond 
to two different inputs for the two-body triplet scattering length. The shaded 
band gives the expected error of a leading order calculation. The existence of 
this correlation was (to our knowledge) first pointed out in Ref. |86j . 

A process relevant to big bang nucleosynthesis is thermal proton capture 
{p+d — >^He-|-7). The calculation of this amplitude is complicated by the presence 
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of Coulomb effects. Yet, a first step towards the goal of calculating such processes 
within the short-range EFT framework was performed in [88ti89j. In this work the 
authors considered thermal neutron capture {n-\-d ^ ^H-l-7) and calculated the 
total cross section at zero energy up to N2L0. The calculation of cross sections 
for this process is easier since no Coloumb effects need to be considered in this 
reaction. The final result is 

(jtot = [0.485(LO) + 0.011(NLO) + 0.007(N2LO)] mb 

= [0.503 ± 0.003] mb , (97) 

where the remaining uncertainty is an estimate of higher order effects. This 
theoretical result is in excellent agreement with the experiment value crf^ = 
0.508 ±0.015 mB. 

A first calculation of Coulomb effects in proton-deuteron scattering (in the 
quartet channel) was done in [90]. However, the energies relevant to big bang 
nucleosynthesis are significantly lower than considered in this publication and 
still a challenge. 



4.3 Four and more Nucleons 

The short-range EFT has also been applied to the four-nucleon sector \91\ . The 
binding energy of the a-particle was calculated using the same approach as dis- 
cussed in subsection l2.4[ In the case of the nuclear system, two effective potentials 
have to be introduced: 

Vs = Xsrs\g){g\ , (98) 
Vt = \tVt\g){g\ , (99) 

where Vs {Vt) denotes the spin-singlet (triplet) potential. Vs and Vt projects on 
the spin-singlet and triplet channels, respectively. The coupling constants As and 
At are therefore renormalized to the singlet and triplet scattering lengths. The 
effective three-body potential is 

V3 = X3Va\0{^\ , (100) 

where Va is the operator that projects on the completely antisymmetrized three- 
nucleon state. The three-body coupling constant A3 can be adjusted using the 
triton binding energy. 

The correlation line between triton and the a particle binding energy were 
calculated using the Faddeev and Faddeev-Yakubovski equations. In Fig. [TSl we 
show this correlation which is also known as the Tjon line. Since the short-range 
EFT uses the minimal number of degrees of freedom for the construction of 
the interaction, it explains therefore the old observation of the linear correla- 
tion between Bt and Ba- In Fig. [15] we show for example also calculations using 
phenomenological potentials [92] and a chiral EFT potential with explicit pions 
[93\ [M] • All these calculations have to lie in the band generated with the short- 
range EFT since these potentials reproduce the large scattering length of the 
internucleon interaction. By the same logic, the short-range EFT also explains 
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Figure 15. The correlation between triton and alpha particle binding energies, as predicted 
by the short-range EFT. The grey circles and triangles show various calculations using phe- 
nomenological potentials [92]. The squares show the results of chiral EFT at NLO for different 
cutoffs while the diamond gives the N^LO result [531 The cross shows the experimental 
point. 



why different versions of the renormahzation group evolved potential l^ow k re- 
produce the Tjon line [97] . 

The short-range EFT was also used to calculate the properties of heavier 
systems than the a-particle. Stetcu et al. [98j considered the a-particle and 
^Li using the short-range EFT in combination with the no-core shell model. 
They reproduced the ground-state results for the a-particle and gave an estimate 
for the binding energy of its first excited 0"*" state. The agreement with the 
experimental value of the ^Li ground state was at the 70 % level. 

4.4 Halo Nuclei 

Another exciting arena for the short-range EFT is the physics of halo nuclei. 
Halo nuclei are weakly bound systems which consist of a tightly bound core and 
a "halo" of one or more nucleons. The radius of a halo nucleus is typically much 
larger than the radius of the core since the total binding energy of such a system 
is comparable to the binding energy of its core. A well-known example is ^He, a 
system consisting of an a-particle core and two neutrons. 

Over the last years ab-initio, wave functions methods (that start from a real- 
istic internucleon interaction) have made considerable progress in describing the 
lighter of the known halo systems [99j , however calculations for heavier systems 
such as ^^Li or ^'^C are not under control yet. The short-range EFT offers an 
alternative perspective on these systems by treating the core of the halo as a 
separate degree of freedom. A system such as ^He becomes therefore an effec- 
tive three-body problem and it's properties can be calculated with the methods 
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Figure 16. Boundary curve in the yjEZ/Bl^ vs. yEnnJE^ plane with leading order error 
bands. Boundary curve shown for a core mass of A = 18 with the experimental data for 
from Ref. [T04l. 



discussed in the previous sections. Such an approach is interesting for several 
reasons. It allows to test whether these systems have features related to large 
scattering length physics such as excited Efimov states but provides also testable 
predictions for observables such matter and charge radii. A further advantage of 
this approach is that scattering observables become directly accessible and can 
be calculated without any further approximations. 

The first application of the short-range EFT to halo nuclei was carried out 
in Refs. |100[ llOlj . In these works the authors considered the one-neutron halo 
^He and calculated in particular phase shifts and cross sections for elastic a- 
nucleon scattering. A further example of a nuclear two-body cluster that has 
been considered is the 2-a system |102j . 

Recently, Canham and Hammer [103j performed the first EFT calculation for 
two-neutron halos. In their work they calculated the binding energies and radii 
of halos such ^^Li and ^''C. Canham and Hammer also addressed the question 
whether any of the considered systems supports an excited Efimov state. Fig. [T6] 
shows a parametric plot {{End B^^'^Y/'^ versus {Enn/ B^'^Y^'^) which describes 
the region in the two-body parameter space that supports a three-body state 
above B^\ They found that the ^''C system might exhibit an excited Efimov 
state close to the threshold. The ^He halo nucleus was not considered by Canham 
and Hammer since it is expected to have P-wave resonance in the neutron-core 
subsystem. This feature raises also the question whether a three-body force is 
required for consistent renormalization in this system. A successful description of 
the effective three-body system would furthermore lead to the exciting problem 
of describing ^He (the effective five-body problem) within the short-range EFT. 
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4.5 Challenges and Opportunities II 

Currently, calculations for observables have only been performed up to N2L0 in 
the EFT expansion. A computation up to N3L0 would provide further informa- 
tion on the convergence radius of the short-range EFT. This is also required to 
fulfill the original promise of outstanding accuracy for few-body observables in 
the few-nucleon sector. 

The inclusion of Coulomb effects has also to be addressed in the near future. It 
is known in the two-body case how to calculate the scattering of, for example, two 
protons. This problem has not been solved to full satisfaction in the three-body 
case although it is crucial for the calculation of scattering processes relevant to 
nuclear astrophysics such pd — > ^He7. Progress in this direction is also required 
for further applications of the short-range EFT in the four-nucleon sector such 
as a calculation of the cross sections for the processes commonly denoted as HEP 
(^He+p ^ ^He+e+ + z^e) and HEN (^He+n ^ '^He+7). 

The calculation of such four-body processes will furthermore require the de- 
velopment of numerical tools which facilitate the calculation of the corresponding 
scattering amplitude^. With such tools at hand, the existence of further univer- 
sal correlation lines could be explored (e.g. the correlation between nucleon-triton 
scattering length and binding energy as presented for example in [106j ). 

5 Final Words 

We have demonstrated that the short-range EFT is an excellent tool to analyze 
the properties of systems whose constituents exhibit a large two-body scattering 
length. 

Atomic physics provides a good testing ground and gives additional justifi- 
cation to analyze the more general (scattering length dependent) implications 
for such systems. It furthermore gives challenges to find solutions to problems 
that are not encountered in the nuclear physics case such as, for example, deeply 
bound two-body states that need to be accounted for in a sound manner. 

The short-range EFT has also provided an alternative perspective on the 
topic of three-body forces. The importance of three-body forces depends strongly 
on the resolution at short distances of the chosen approach. In EFT approaches 
to nuclear systems that employ pionic degrees of freedom, the first three-body 
forces are generally required at higher orders in the EFT expansion. There is 
therefore nothing fundamental about a three-body force, but it is rather a further 
necessary tool to account for unknown short-distance physics. 

It has also been shown that pionic degrees of freedom are irrelevant to the 
description of a number of important quantities such as the triton bound state 
or the a-particle. A number of observables have been computed beyond leading 
order and excellent agreement with experiment and/or calculations using realistic 
inter nucleon was found. 

We have tried to supply a number of open problems that should be addressed 
in the near future. We hope they are understood not as barriers but as excep- 



^While this manuscript was finalized we became aware of a first calculation of four-nucleon 
scattering observables using the framework of the short-range EFT [105] . 
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tional chances for the short-range EFT to largely extend its current region of 
applicability. 
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